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The magnetization of Pbi_a:Eua;Te samples with x = 1.9, 2.6 and 6.0% was measured at 20 mK 
in fields up to 50 kOe, and at 0.6 K in fields up to 180 kOe. The 20 mK data show the magnetization 
steps (MSTs) arising from pairs and from triplets. The pair MSTs are used to obtain the dominant 
Eu-Eu antiferromagnetic exchange constant, J/ks = —0.264 ± 0.018 K. The exchange constant for 
triplets is the same. Comparison of the magnetization curves with theoretical simulations indicates 
that the Eu ions are not randomly distributed over all the cation sites. The deviation from a random 
distribution is much smaller if J is assumed to be the nearest-neighbor exchange constant Ji rather 
than the next-nearest-neighbor exchange constant J2. On this basis, J is tentatively identified as Ji. 
However, the possibility that J = J2 cannot be excluded completely. To obtain agreement with the 
data, it must be assumed that the Eu ions tend to bunch together. Comparision with microprobe 
data indicates that the length scale for these concentration variations is smaller than a few /im. The 
theoretical simulations in the present work improve on those performed earlier by including clusters 
larger than three spins. 
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I. INTRODUCTION 

An important group of dilute magnetic semiconductors 
(DMS) are lead salts in whicL a fraction of the Pb ions 
have been replaced by Eu ions.u Pbi_a:Eua:Te is one mem- 
ber of this group. The other members are Pbi_a:Eua;Se 
and Pbi_2:Eua;S. All these materials have the rock-salt 
structure, with an fee cation lattice. The Eu^+ ion, with 
7 electrons in the half-filled 4/ shell, has zero orbital an- 
gular momentum and total spin S=7/2. EPR data show 
that in Pbi_-i;Eua;Tc g = 1.98 and that the spin hamilto- 
nian for the Eu^+ ion contains a very small crystal-field 
anisotropy.B 

It has been known for several years that the Eu-Eu 
exchange interaction in these IV- VI DMS is two orders 
of magnitude smaller than the Mn-Mn exchange interac- 
tion in the traditional II- VI DMS (e.g., Cdi-^^Mn^^Te). 
In both types of DMS the exchange interaction is an- 
tiferromagnetic (AF), but the leading Eu-Eu exchange 
constant J is of order —0.1 K, compared to —10 K for 
the Mn-Mn exchange constant. 

Much of the early information concerning these Eu-Eu 
exchange interactions came from measurements of the 
Curie- Weiss temperature 9, and from ajialysis of high- 
field magnetization data at T = 4.2 K.tTEl As discussed 
later, both of these methods yield only a rough esti- 
mate of the dominant AF exchange constant J. In the 
present paper we present an accurate determination of 
J in Pb]pjrjEua;Te using the magnetization-steps (MSTs) 
method.EHO Because \ J\ is of order 0.1 K, the observation 
of the MSTs required the use of a dilution refrigerator 
operating well below 0.1 K. 

In addition to yielding J, the low-temperature magne- 
tization data also gave information about the distribution 
of the Eu ions over the cation sites . Such information is 
not readily available by other means. The analysis indi- 
cates that the distribution is not perfectly random (equal 
probability of occupation of all cation sites in the crystal) . 
Instead, the Eu ions tend to bunch together. Another 
important issue is whether the dominant AF exchange 
constant J corresponds to the nearest-neighbor (NN) or 
to the next-nearest-neighbor (NNN) exchange constant, 
Ji or J2 , respectively. This issue is addressed on the 
basis of the results for the Eu spatial distribution. 

II. THEORY 

Much of the relevant theory was discussed in a recent 
paper on the MSTs in Pbi_a;Eua;Se.El We refer the reader 
to this earlier paper, and confine ourselves here to a sum- 
mary of the main points. Some newer theoretical results 
are also mentioned. 

As a first approximation we consider a simple model 
which captures the main features of the experimental 
data. In this "single- J" model for a DMS, only one AF 



exchange constant J is included. Other exchange con- 
stants, and all anisotropies, are neglected. In this model 
each Eu ion on the fee sublattice belongs to a particular 
type of "cluster". The cluster types considered explic- 
itly in Ref. ^ are: singles (isolated magnetic ions which 
are not connected by any exchange bonds), pairs, open 
triplets (OTs), and closed triplets (CTs). More recently 
the sipc types of quartets in the fee lattice were considered 
also.0 

Assuming that J is negative (antiferromagnetic) , the 
magnetization curve at low T {ksT « \J\) for clusters 
of a given type has two main features: a quick saturation, 
of the zero- field moment followed by a series of MST's.&lI 
This behavior is illustrated by the exact calculation of 
the magnetization curve for small clusters in this simple 
model. We will use normalized fields /i„ = g^BHn/\J\ for 
the step positions. Pairs of Eu^+ ions {S — 7/2) have no 
zero-field moment. The seven MSTs from pairs occur at 
normalized fields /i„ = 2, 4, ...14. Thus the pairs saturate 
at qhbH = 14| J|. The next most important clusters are 
open triplets. They have a zero-field groundstate with 
a net spin of 7/2 which saturates quickly. The MSTs 
occur at /i„ — 9, 11. ..21. Closed triplets have their MSTs 
at hn = 1,3. ..21. Note that the last three MSTs due to 
triplets occur when the pairs are already saturated. For 
five of the six types of quartets the series of MSTs end 
at 28| J|; only the "stung quartet" ramp ends at a lower 
field gfigH = 24.2| Jj.E 

The exchange constant J is usually obtained from the 
observed Hn for pairs since, except for singles, these are 
the most numerous. As was done in Ref. |^, we esti- 
mate an uncertainty in the determination of J due to 
simplifications introduced by the single-J model. The 
cubic crystal-field anisotropy, the dipole-dipole interac- 
tion, and exchange constants other than J basically shift 
and broaden the MSTs. Their effects on the positions of 
the MSTs from pairs are included in the equation 

gfiBHn = 2n\J\ + Ar, (1) 

with n = 1,2. ..7. The shifts A„ depend on n, the orien- 
tation of the sample (anisotropy) and the concentration 
(further neighbor exchange constants). We have calcu- 
lated the magnetization and step postilion of pairs taking 
into account the cubic crystal field arid dipole-dipole in- 
teractions, using known parameters.l3 If J is determined 
from a fit to Eq. ([|), assuming a constant A„ (the pro- 
cedure used below) , the resulting error due to the crystal 
field anisotropy is less than 5%. The error due to dipole- 
dipole interactions is less than 1%. 

Exchange interactions with further neighbors (not in- 
cluded in tiie single-J model will also shift and broaden 
the MSTs.t2l These effects will be discussed later m con- 
nection with the data analysis. 

The magnetization M of the sample as a whole is ob- 
tained by adding the contributions of the various cluster 
types. The contribution of each cluster type is the prod- 
uct of the magnetization per cluster and the number of 
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clusters of that type. To calculate the number of clusters 
of a given type one needs to know how the magnetic ions 
are distributed over the cation sites. Normally, a random 
distribution is assumed. The popwlations of the various 
cluster types are then well knownJj'til As discussed later, 
deviations from a random distribution can be detected 
by analyzing the measured magnetization curve at low 
temperatures. The ability to detect such deviations is a 
major advantage. Any tendency of the magnetic ions to 
bunch together or to avoid each other is revealed. The 
determination of the AF exchange constant J is indepen- 
dent of the spatial distribution however. 

In Ref. ^ the theoretical simulations of the magnetiza- 
tion curves included clusters up to triplets. In the present 
work the simulations were improvedJpy adding the contri- 
butions of the six types of quartets. □ In addition a rough 
estimate of the contribution of clusters larger than quar- 
tets was also included, as will be discussed shortly, so 
that all the spins were accounted for. 

Clusters larger than quartets (quintets, sextets, etc.) 
will be referred to collectively as "others" . If the zero- 
field ground state of any such cluster has a net spin then 
this net spin will align rapidly with H at low T. At 
higher fields a series of MSTs will occur. In practice, the 
MSTs from large clusters are very small in size, so that 
they are not resolved. The series of MSTs from a given 
cluster type then merges to form a ramp. The ramp ends 
at a field which depends on the cluster type. 

The magnetization curve of the others is therefore a 
sum of many initial fast rises of M followed by a super- 
position of many ramps ending at different fields. Here, 
we use an approximation in which the contribution of the 
others is represented by a single initial fast rise of M fol- 
lowed by one ramp. The initial fast rise is approximated 
by a Brillouin function for spin 7/2 with a saturation 
value corresponding to 1/5 of the saturation value of the 
others. The 1/5 weight is motivated by earlier results for 
the initial rise of the magnetization.El The remaining mag- 
netization rise (4/5 of the saturation value of the others) 
is approximated by a single ramp which starts at _ff = 
and ends at g^sH — 35| J|. The latter value is expected 
to be slightly higher than the average saturation value for 
all quintets but lower than the average saturation value 
of sextets. 

For our samples, the contribution of others to the mag- 
netization was small, only 0.2% of the saturation mag- 
netization Mq for X = 1.9%, and only 0.6% of Mq for 
X = 2.6%. Even for the sample with x — 6% the others 
contributed only 9% of Mq. Our approximation for the 
contribution of the others should therefore be adequate. 

III. EXPERIMENTAL TECHNIQUES 

The three Pbi_a;Eua;Te samples were grown by the 
Bridgman method. The Eu concentration x was deter- 
mined from the saturation magnetization. A moment of 



6.93^_B per Eu ion (based on S* = 7/2 and g = 1.98) was 
assumed. A small correction for the lattice susceptibility, 
Xd — Sx 10^ emu/g, was applied.0 The results for the 
three samples gave x — 1.9, 2.6 and 6.0%. These val- 
ues are supported by the Curie constants, obtained from 
susceptibility data, which gave x = 1.9, 2.4, and 5.9%, 
respectively. We shall adopt the first set of values, which 
we regard as more accurate. 

The Eu concentration x was also obtained from micro- 
probe measurements. Approximately 35 spots on a single 
surface of each of the samples were probed. Each spot 
had a diameter of 5 /im, and probed the Eu concentration 
to a depth of about 2 /zm. The average values for x and 
the standard deviations were 1.7 ± 0.3%, 2.5 ± 1.1%, and 
5.6 ± 0.2%. The large standard deviation for the second 
sample is due to a small region of higher concentration 
near one corner. The other two samples, particularly the 
one with the highest concentration, were quite homoge- 
neous on length scales larger than several /im. 

Magnetization measurements at 20 mK were made us- 
ing a force magnetometer which operated in the mixing 
chamber of a plastic dilution-, pefcigerator. This equip- 
ment was described earlier .I3l9't2lll3 The main magnetic 
field H, up to 50 kOe, was generated by a superconduct- 
ing magnet. The force was produced by a superimposed 
dc field gradient, dh/dz = 0.8 kOe/cm, which was gen- 
erated by an independent set of superconducting coils. 
With this gradient the variation of the field over the vol- 
ume of the sample was less than 0.2 kOe. None of the 
samples were oriented, so that the direction of the field 
relative to the crystallographic axes was not known. As 
discussed later, the anisotropy in Pbi-^EuajTe is small, 
so that field orientation is not critical. 

The magnetization was also measured in fields up to 
180 kOe using a vibrating sample magnetometer (VSM) 
which was adapted for use in a Bitter magnet. The sam- 
ples were immersed in liquid ^He at 0.6 K. Other mag- 
netization data were taken with a SQUID magnetometer 
systemtJ at temperatures above 2 K and in fields up to 
50 kOe. Among the data taken with this system were 
the low-field susceptibility data used to obtain the Curie 
constant and the Curie- Weiss temperature 6. 

IV. RESULTS AND DISCUSSION 

A. Magnetization curves 

Figure |l| shows the top half of the magnetization curve 
at 20 mK for x = 1.9%. The fast rise at the lowest fields is 
due to the alignment of singles, with minor contributions 
from other clusters with zero-field moments. This initial 
fast rise of M is followed by a large ramp on which six 
MSTs are clearly visible. The six MSTs are due to pairs. 
Only six, instead of seven, MSTs are visible because the 
first MST is masked by the large initial fast rise of M on 
which it is superimposed. The ramp due to pairs ends 
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FIG. 1. Top half of the magnetization curve at 20 mK for 
the sample with x = 1.9%. The magnetization M has been 
normalized to its saturation value Mo. 

near 28 kOe. It is followed by a much less steep ramp 
due to open triplets (OTs) which ends near 43 kOe. (The 
OT ramp is predicted to start well before the pair ramp 
ends, but because the pair ramp is so much steeper the 
OT ramp does not stand out in this field region.) For this 
sample, with the lowest concentration, the magnetization 
becomes practically saturated once the OT ramp ends. 

The lowest curve in Fig. § shows the differential suscep- 
tibility (field derivative of the magnetization) obtained 
numerically from the data in Fig. |l|. The six MSTs due 
to pairs appear as large peaks. Three much smaller peaks 
are seen between 34 and 43 kOe. These are the last MSTs 
from the OTs. To our knowledge this is the first clear ob- 
servation of MSTs from triplets in a DMS. The other two 
curves in Fig. ^ show similar results for pair and triplet 
MSTs in the other two samples. 

Figure ^ shows magnetization curves for two of the 
samples at 0.6 K. These curves extend to fields of about 
180 kOe. Also shown, for comparison, are the corre- 
sponding magnetization curves at 20 mK, adjusted so 
that they agree with the 0.6 K data at 50 kOe. The main 
differences between the 0.6 K and 20 mK data are that at 
0.6 K the initial rise of M is more gradual, the pair ramp 
is more rounded, and the pair MSTs are not resolved. 
These differences are the expected effect of temperature. 
The data at 0.6 K show that for x — 1.9% complete sat- 
uration is achieved near 60 kOe. For x = 6.0% however, 
the approach to saturation is more gradual because a 
larger percentage of the spins are in large clusters which 
saturate more slowly. 
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FIG. 2. Field derivative of the normalized magnetization 
m — M/Mo for the three samples at 20 mK. These results 
were obtained by a numerical differentiation of the magneti- 
zation curves. 

B. Dominant AF exchange constant J 

In all three samples the pair and triplets involving the 
exchange constant J gave rise to distinct ramps on which 
well-resolved MSTs were observed. No additional ramps 
or MSTs from any other exchange constant were found. 
The single- J model therefore seems to be appropriate for 
describing the magnetization curves. Additional evidence 
that the corrections for the single- J model are small will 
be discussed later. 

The fields iJ„ at the pair MSTs were obtained from 
the peaks in Fig. ^. Figure |^ shows a plot of these fields 
for the three samples. Because the first MST was not 
observed, the plot is only for n = 2-7. It is evident that 
Hn is approximately linear in n, so that Eq. ( p]) w ith a 
constant A„ is a good approximation (see Ref. |15|). 

The values for J obtained from least squares fits to 
Eq. (0), holding A„ constant, were the same for all three 
samples within 2%. The average value was J/kB = 
—0.264 K. The values for A„ obtained from fits of the 
three data sets to Eq. (^ were all positive. The aver- 
age was A„ = 0.69 kOe, which is only a fraction of the 
spacing AH — 4 kOe between successive MSTs. 

As far as further neigbor exchange interactions are con- 
cerned, we can say the following. The data show that 
any other AF exchange constant is at least a factor of 3 
smaller than the observed J. If this were not the case, an- 
other noticeable ramp ending between 10 kOe and 30 kOe 
should have been present, leading to a change of slope 
in this field region. Such a change of slope was not 
observed. Furthermore, the observed MSTs show only 
a small amount of non-thermal broadening and A„ de- 
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FIG. 3. Magnetization curves for x = 1.9 and 6.0% at 
0.6 K. These results have been corrected for the lattice dia- 
magnetism. Also shown are the magnetization curves at 
20 mK, with the vertical scale adjusted so that they match 
the 0.6 K data at 50 kOe. 



pends only weakly on x. These observations show that 
the magnetization above about 10 kOe is very little af- 
fected by other (smaller) exchange constants originating 
from further neigbors. 

The following arguments indicate that the quoted value 
J/ks — —0.264 K corresponds to the largest AF ex- 
change constant. The 0.6 K data in Fig. ^jshow no other 
MSTs or ramps in fields up to 175 kOe. This means that 
if there were any larger AF exchange constant, its magni- 
tude should have exceeded 12 K. On the other hand the 
Curie- Weiss temperatures discussed below rule out such 
a large AF exchange constant. Thus, the observed MSTs 
gave the largest J. 

Our final result is then that the dominant exchange 
constant has the value J/ke = -0.264±0.018 K. The 7% 
uncertainty is motivated by the calculations of the cor- 
rections mentioned in Sec. and the observation that 
further neighbor interactions affect the MSTs very lit- 
tle. This J for Fbi-^^EuajTe is only slightly larger than 
—0.24 K for Pbi_2:Eu2:Se.cl There was no measurable de- 
pendence of J on a; in the present Pbi_a:Eua;Te samples 
(see Fig. m, despite the fact that the band gap changes 
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FIG. 4. The fields Hn at the MSTs from pairs as a func- 
tion of the step number n. The straight line represents Eq. (|^) 
with the average values for J and A„. 

by nearly a factor of 2 as a; changes from 1.9 to 6.O%.0 
The fields at the three observed MSTs from OTs, above 
32 kOe, are close to those expected from the value of J de- 
rived from the pair steps. This means that the exchange 
constants for pairs and triplets are the same (within a 
few percents), as expected. 

Low-field susceptibility data were taken at tempera- 
tures down to 2 K. The most accurate result for 9 was fpt 
the sample with x = 6.0%. With the usual assumptionsll3 
the value 9 = —1.9 K gave J/fcs = —0.25 K (assuming 
J = Ji). The results for the other two samples, with 
lower X, were J/ks = -0.30 {x = 2.6%) and -0.22 K 
{x = 1.9%), but these were judged to be less accurate. 
The values found by Gorska et aJ.cl were J/ks — —0.38 K 
for X = 3%, and —0.27 K for x = 6%. The spread in the 
values of J obtained from 9 shows that that the exper- 
imental uncertainty is not negligible, especially for low 
concentrations. Determining J from the Curie- Weiss 9 
has several other drawbacks. First, it assumes that the 
distribution of the Eu ions is random, which, as discussed 
later, may not be exactly true for Pbi_2,Eu2,Te. Second, 
9 depends on all exchange constants. Third, even when 
one J is much larger than all the others, its estimate de- 
pends on whether it is identified as Ji (between NNs, as 
was done above) or J2 (between NNNs). The identifica- 
tion of J as J2 would lead to an estimate of J which is 
larger by a factor of 2 compared to J = Ji . The MSTs 
method is a direct determination of J, which is indepen- 
dent of the identification of J, and is also independent of 
the spatial distribution of the Eu ions. 

Gorska et ala also determined J by analyzing the mag- 
netization curve at 4.2 K. At this relatively high temper- 
ature, fc^r = 16|J|, the MSTs are not resolved, and the 
ramps due to pairs and OTs do not stand out clearly. 
The analysis which was performed fitted the magnetiza- 
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FIG. 5. Comparison between the measured magnetization 
curve at 20 mK with computer simulations based on the Ji - 
and J2 models, a) a; = 1.9%. b) a; = 2.6%. 
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FIG. 6. Comparison between the measured magnetization 
curve for x — 6.0% and computer simulations based on the 
Ji - and J2 models, a) T = 20 mK, b) T = 0.6 K. 



tion curve to a sum of two contributions: one from sin- 
gles and the other from "PAIRS" . The value for J was 
deduced from the PAIR contribution. The results were 
J/kB = -0.43 K for X = 3%, and -0.50 K for a; = 6%. 
The reason that these values are much higher than our 
value of —0.264 K, is that the assumed "PAIRS" actually 
included not only true pairs but also triplets, quartets, 
and larger clusters. Since the saturation field of a cluster 
increases with cluster size, the net effect was that the as- 
sumed PAIRS saturated at higher fields than true pairs. 
This caused the deduced J to be higher than the true J. 
It is probably significant that the fitted J was larger for 
X — 6% than for x = 3%, because the percentage of large 
clusters increases with x. 



C. Eu distribution and the identity of J 

The identity of the dominant AF exchange constant 
J, whether it is Ji for NNs or J2 for NNNs, is a sig- 
nificant issue. Normally, exchange constants decrease 
rapidly with distance, so that the dominant J is between 
NNs. This expected noririal behavior has been found in 
all Mn-based II-VI DMS.H On the other hand, in pure 
EuTe the largest exchange constant is J2, which is an- 
tiferromagnetic with a value of about —0.2 K. The ex- 



change constanLpL in EuTe is ferromagnetic and smaller 
in magnitude.EZrEa Thus, after J was measured in the 
present work, it was still not obvious whether it was Ji 
or J2. 

To address this problem the measured magnetization 
curves were compared with computer simulations based 
on the two competing hypotheses: J — Ji (Ji model), 
or J = J2 (J2 model). Both models are single- J models. 
The simulations were similar to those in Ref. || except 
that they also included the quartets and the others. We 
also chose to normalize the magnetization curves to the 
saturation value Mq rather than to the data point at the 
highest field (see Ref. |9|). The key assumption in the 
simulations was that the Eu ions were distributed ran- 
domly, i.e., the probability of occupation of each cation 
site in the crystal was x. The simulations with the Ji 
- and J2 models lead to different magnetization curves, 
essentially because there are 12 NNs but only 6 NNNs 
in the fee cation lattice. This difference means that in 
the J2 model the number of singles is higher and the 
populations of pairs and larger clusters are lower. 

Figure || compares the magnetization curves at 20 mK, 
for X — 1.9 and 2.6%, with simulations based on the 
two models. To account for the somewhat larger ob- 
served width of the MSTs than that expected from ther- 
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mal broadening alone, the simulations of data at 20 mK 
used an effective temperature of 100, mK. This change 
has no effect on the discussion below .B The simulation of 
data at 0.6 K were made with the actual temperature. 
From the figure it is clear that the Ji model simulation 
is closer to the data, but even this model fails to repro- 
duce the measured curves exactly. The deviations are 
more pronounced for x = 2.6%. The observed initial 
rise of M is smaller than predicted, indicating that even 
the Ji model overestimates the number of singles. Fur- 
thermore, the observed slope between 30 and 40 kOe is 
higher than predicted by both models. This means that 
there are more triplets and/or larger clusters than calcu- 
lated using a random distribution. The same features are 
also seen in Fig. which is for x = 6.0% at 20 mK and 
at 0.6 K. The Ji model is closer to the data, but even 
compared to this model the actual number of singles is 
lower and the number of triplets and/or larger clusters is 
higher. 

These results indicate that the local Eu concentration 
in the vicinity of a typical Eu ion is higher than the av- 
erage concentration for the sample as a whole. Thus, the 
Eu ions tend to bunch together, leading to an inhomoge- 
neous Eu distribution. The size of the regions with higher 
Eu concentration cannot be inferred from analysis of the 
magnetization curves alone. The microprobe measure- 
ments (Sec. [II) examined the concentration variations 
on length scales greater than a few fun. The attempt to 
account for the magnetization curves by using the mi- 
croprobe results as the local concentration profile was 
unsuccessful. This is particularly evident for the sample 
with X = 6.0%, in which the microprobe concentration 
showed very little variation with position. Combining the 
analysis of the magnetization data with the microprobe 
results we conclude that the length scale of the concen- 
tration inhomogeneities was smaller than a few /im. 

A quantitative measure of the degree to which the Eu 
ions bunch together can be obtained by introducing the 
concept of a local concentration xl. The simple picture 
which is behind this concept presumes that there is a 
typical local concentration in the vicinity of a typical Eu 
ion. It can differ from the average concentration x for 
the sample as a whole. The populations of the various 
clusters are calculated assuming that the probability of 
occupation of each cation sites in the vicinity of a Eu ion 
is XL instead of x. The value of xl is found by varying 
the Eu concentration in the simulation until a match with 
the experimental magnetization curve is achieved. Using 
this procedure the following results were obtained for the 
samples with x = 1.9, 2.6 and 6.0%. With the Ji model, 
XL = 2.5, 3.8, and 7.7%, respectively. With the J2 model, 
XL — 4.7, 7, and 14%, respectively. Obviously, a much 
smaller difference between xl and x is required to achieve 
agreement with the Ji model. 

In the discussion above, we have assumed that the de- 
viations of the data from a single-J model were due to 
the assumption of a random distribution. An alternative 
would be to question the validity of the single-J model 



itself. It is known that long range exchange interactions 
between further neighbors can influence the maEp|etiZ|a-| 
tion process in diluted magnetic semiconductors. t3'E2rEHl 
Such long-range interactions may be important in small- 
gap DMS, such as the present system. The key ques- 
tion here is whether long-range interactions could have 
affected the magnetization in fields above 10 kOe signifi- 
cantly, since the conclusion of a non-random distribution 
was based on the data at these high fields, especially 
above 30 kOe. 

The effects of further neighbors (distant neighbors) on 
the magnetization curve have been treated using sev- 
eral approximate methods. One approach starts from 
the clusters in the single-J model, but then subjects 
these clu^rs to effective fields arising from further 
neighbors.tj Alternative approaches use more general 
spin clusters which include further neighbors within the 
clusters, so that the intracluatec—interactions already 
include the further neighbors.E3L3 Since our analysis 
started from the single-J model, the first approach (ef- 
fective fields acting on the single-J clusters) is more con- 
venient here. The effective fields from further neighbors 
slow down the alignment of spins which are singles, and 
they shift and broaden the ramps and MSTs arising from 
larger clusters (e.g. pairs and triplets). These effects in- 
crease with the Eu concentration x. We now present sev- 
eral arguments which indicate that in the present work 
the further-neighbor effective fields were far too weak to 
significantly affect the magnetization curve well above 10 
kOe. 

At 20 mK the magnetization rises very quicly at low H 
(Figs. |l] and ||). The singles seem to become practically 
saturated at about 2 kOe. Considering that the cubic 
crystal field anisotropy also slows down the alignment 
of the singles ,eI we conclude that typical further- neighbor 
effective fields acting on singles are less than about 1 kOe. 
Such weak effective fields are unlikely to have a significant 
effect on the magnetization well above 10 kOe. 

The fields Hn at the MSTs from pairs are shifted by 
the further-neighbors effective fields.ll2l These shifts, in- 
cluded in the A„ of Eq. (|l|), are predicted to increase 
with X. However, the results in Fig. ^ indicate that the 
A„ are all less than about 1 kOe. Also, despite the 
factor-of-three change in x, the fields Hn do not vary 
by more than a fraction of 1 kOe. This means that typ- 
ical further- neighbor effective fields acting on pairs were 
no more than a fraction of a kOe, even for the sample 
with X = 6%. Typical effective fields acting on triplets 
and quartets should not be higher by more than a factor 
of 2 or so. It is very unlikely that such weak effective 
fields will have a significant effect on the magnetization 
curve well above 10 kOe. 

The spread in the magnitude of the effective field acting 
on different-i)airs gives rise to a broadening of the MSTs 
from pairs. t3 Such a broadening should increase with x. 
The fact that well resolved MSTs were observed even for 
X — 6% means that the spread in the effective fields was 
considerably smaller than the 4 kOe separation between 
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the MSTs. Thus not only was the average effective field 
small, but the spread was also small. Such a distribution 
of effective fields should not have a significant effect on 
the magnetization curve well above 10 kOe. 

The final argument in support of a non-random distri- 
bution is based on an analysis of the change in the slope 
of the magnetization curve near 43 kOe. This change oc- 
curs when the triplet ramp ends. (The change in slope is 
visible in Figs. |^ and||, but is clearer in expanded plots of 
the magnetization curves at 20 mK.) The magnitude of 
this change in slope is related to the number of triplets. 
Analysis of the data in all three samples indicates a sig- 
nificantly larger number of triplets than predicted assum- 
ing a random distribution. Further-neighbor interactions 
cannot account for the significanly larger change in slope 
which was observed. 

Claims of a non-random distribution in 11- VI DMS 
have been made more than a decade ago (see e.g. 
Ref. ^). However, MST studies showecLthat the distri- 
bution in such DMS was in fact random.cJ'El Here, in our 
Pbi_a;Eua;Te samples, we find a small non-randomicity 
which we believe to be genuine. In Ref. |^ (Pbi_^Eu^Se) 
deviations of the data with the Ji-model were also found, 
but the simulations included only clusters up to triplets. 
A re-analysis of those data with our improved model, 
which includes quartets and an approximation for larger 
clusters, still does not lead to perfect agreement. As 
in the present work, good agreement with the data is 
obtained if small deviations from a random distribution 
are assumed (and that J = Ji). The non-randomicity 
in the Pbi_2:Eua;Se samples is smaller than that in the 
Pbi_a;Eua;Te samples. 

Returning to the issue of the identity of J , there are 
two possibilities. Either J is Ji in which case the Eu ions 
have only a fairly modest tendency to bunch together, or 
J is J2 with a very large increase of the local Eu concen- 
tration. Because the jusual distribution of magnetic ions 
in a DMS is random,El we believe that the first possibil- 
ity is more likely. On this basis we tentatively identify 
J = —0.264 ± 0.018 K as Ji. However, since the actual 
distribution is unknown, the possibility that it is J2 can- 
not be excluded entirely. If one accepts that J is Ji then 
the exchange constants in Pbi-^^Eu^Xe when x is low are 
very different from those in EuTcEZliia Such a difference 
may be the result of a different band structure and a 
different position of the Eu levels. 
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